Intermediate wave-function statistics 
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We calculate statistical properties of the eigenfunctions of two quantum systems that exhibit 
intermediate spectral statistics: star graphs and Seba billiards. First, we show that these eigen- 
functions are not quantum ergodic, and calculate the corresponding limit distribution. Second, we 
find that they can be strongly scarred, in the case of star graphs by short (unstable) periodic orbits 
and in the case of Seba billiards by certain families of orbits. We construct sequences of states which 
have such a limit. Our results are illustrated by numerical computations. 



It has been conjectured that the quantum spectral 
statistics of systems that are chaotic in the semi-classical 
limit are generically those of Random Matrix Theory Q . 
The behaviour of the eigenfunctions of such systems is 
described by the semi-classical eigenfunction hypothesis 
0, 1^, which implies that they equidistribute over the 
appropriate energy shell. This is in agreement with a 
theorem of Schnirelman 4] which implies equidistribu- 
tion of almost all eigenstates on scales independent of 
h, assuming only classical ergodicity. Such behaviour is 
termed quantum ergodicity. This theorem still permits 
the possibility of a small number of states which do not 
equidistribute. 

It has been suggested that some of these exceptional 
states may be "scarred" by short classical periodic orbits 
Further investigations d H H E have distin- 
guished between weak and strong scarring. Weak scar- 
ring relates to states averaged over energy windows that 
contain a semi-classically increasing number of levels, 
whereas strong scarring means that sequences of states 
can be constructed whose limit is wholly or in part sup- 
ported by one-or-more periodic orbits. So far the only 
systems known rigorously to support strong scarring are 
the cat maps 12\ , which have non-generic spectral statis- 
tics dii. 

For systems that are classically integrable, it is ex- 
pected that the quantum spectral statistics are Poisso- 
nian, i.e. those of independent random numbers [l^ . 
The corresponding eigenfunctions semiclassically equi- 
distribute on tori in phase space [isj . 

Recently, classes of systems which exhibit spectral 
statistics that are intermediate between Random Matrix 
and Poissonian have been discovered 0, 0, 0| . Two 
representative families of examples are Seba billiards 
and star graphs [i^l . It was shown in jij that these two 
systems have the same (intermediate) spectral statistics. 
We study the eigenfunction statistics of such systems. 
Specifically, given that these systems are not classically 
ergodic, we are interested in whether the eigenfunctions 
are quantum ergodic, and whether they show strong scar- 
ring (that they exhibit weak scarring may be shown using 
the methods of 



Star graphs are quantum graphs [22| which have one 
central vertex, and 6 outlying vertices each connected 
only to the central vertex [20|. For such graphs, the limit 
& — > oo is analogous to the semi-classical limit. To inves- 
tigate the possibility of quantum ergodicity in this limit, 
we consider a graph with b = av bonds, where w 3> 1, 
a > 1, and introduce the observable B defined by 



B 



1 on bonds indexed 

on bonds indexed u -f 1, . . 



Thus B picks out a fraction a^^ of the bonds. Let i/'n 
denote the wave-function associated with the n**^ eigen- 
state. We calculate the probability distribution, P{R), 
for n chosen at random, that (V'n|i?|V'n) is less than R, 
subject to some mild restrictions on the bond lengths. A 
system which exhibits quantum ergodicity would have 



P{R) 



0, 
1, 



< i? < 
a-^ < R< 1. 



Our result (see equation 10) and figure ^ below) differs 
from this, proving that star graphs are not quantum er- 
godic. In fact we are able to say more: for a fixed (fi- 
nite) number of bonds, we explicitly find eigenstates that 
are strongly scarred along closed (unstable) orbits of the 
graph with period 2. This is the first class of examples 
showing generic (in this case intermediate) behaviour in 
which strong scarring has been rigorously demonstrated. 

The term Seba billiard refers to any integrable quantum 
system that has been perturbed by the addition of a point 
singularity. We consider the specific example of a billiard 
on a torus. By exploiting the connection between Seba 
billiards and star graphs 21] we argue that Seba billiards 
are also not quantum ergodic and find states that appear 
to show behaviour analogous to strong scarring, in this 
case by families of orbits. 

We begin by describing how to calculate the probabil- 
ity distribution P{R). 

Eigenenergies of a star graph with b bonds are given 
by En — where fc„ is the n*'' solution of Z{k) — 
with 



Z{k) = ^tanfcLj 



(1) 



2 



the individual bond lengths being denoted hy Li, Lb. 
The component of the n*^ wave-function on the i^^ bond 
of the graph is ipn,i{x) — Ai(fc„) cos fc„(x — Li), where 



72 



cos knLi{J2j Lj sec2 knLjY/'^ ' 
the sum being taken over all bonds. Then 

sec^ knL 



(2) 



+ 0(fc-i). (3) 



To calculate the distribution of values taken by this quan- 
tity we average over a large number of states, making the 
error term in Q negligible. We choose incommensurate 
bond lengths from an interval [L, L + AL] that shrinks 
in such a way that vAL as w ^ oo. Thus we can 
replace Li by L wherever it does not multiply fc„. 

To evaluate a function f{k) at the zeros of Z{k) we 
integrate against the density of states, so 



1 

N 



N 



n=l 



N 



f{k)Z'{k)5[Z{k)]dk 



where 5 denotes the Dirac delta function. Writing 
the delta function in Fourier representation, 5{x) — 
(27r)-i e'^^dC, and taking the Hmit TV ^ oo. 



1 ^ 

n— 1 

= lim — 

2T:d K^oo A 



(4) 



f{k)Z'{k) exp[iCZ(fc)]dCdfc, 



writing K — k^ and using fcjv ~ N/d, where d = bL/n 
is the mean density of states. We apply l@J with f{k) = 
exp(i/3X^(fc)) where 



X^{k) 



b 



-77 / sec^ kLj 

J=t)+1 



4E 



sec^ kLi (5) 



for /?, J] constants. This is related to the distribution of 
('0„ I i? !?/'„) by the fact that 

P(X^(fc„) > 0) = P((V^„|S|i/;„) < R) 

when i? and r\ are related by?7 = l/i?— 1. 

We observe that k only appears in multiplied by a 
bond length, and as an argument of a 7r-periodic func- 
tion. Since the bond lengths are incommensurate, the k 
integral can be re-written as a multiple integral over the 
6 variables Xj = kLj. A similar argument was used in 
[2^13. The integrand then factorises, so that 



where 



1 



i/377 



/i = — / sec a;exp — tana; 5- sec x dx, 



,'iC, 2 1 , 

exp [ — tan a; — sec x ) dx, 



I3 is obtained by replacing (3 with —P/rj in /2, and /4 
by making the same substitution in /i. Techniques to 
analyse the asymptotics of these integrals were discussed 
in ^] . Using them we find that 



h 



and 
I2 ~ exp 



-^v/i^exp(g)-Cerf(^ 



as u — > cxD. Substituting the above into ((HJ and denoting 
the result e(/3), we arrive at 



1 



1 



exp 



VP 



dC, 



where 
and 

^(0 - - 



1 fie 

exp 



2 ^ fie 



■ exp 
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^ erf 



2y/ir] 



exp 



-^j+C(a-l)erf ^ 



iir/4^ 



The Fourier transform of e(/3) is the probability den- 
sity function of X^(fc„) where the index of the state, n, 
is chosen at random. The probability distribution for 
{ilJn\B\4'n) to be less than R is then given by 



P{R) = — / e{P)e-''^''dpda 
27r Jo J -00 

The Fourier transform of e(/3) is 



(7) 



?)=i/fl-i 



27ra 



-5R 



2\/^ 



having made the substitution ^ = C/v^ using the no- 
tation k;(z) = erfc(— iz). Performing the cr-integral 
in 10 gives, finally, 



lim 



1 ^ 1 /• 

n— 1 



7- riJ— 1 tQ-U — U 

-'1^2 



+ (a - l)hP,ir—'dC (6) 



P(i?) 



T(c)iog(T(e))dc 



(8) 



with r/ = - 1, for < i? < 1. 
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Numerical Study 
"I - Analytical prediction ■ 
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FIG. 1: Comparing P{R), as given by to a direct numer- 
ical computation for a star graph with 90 bonds when a = 3. 
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FIG. 2: Difference between P[R) and numerics when h 
15(+), 30(x), 45(*), 60(a), 75(B), 90(0). 



The results of numerical computations which support 
this calculation are shown in figures ^ and [21 

We now turn to constructing sequences of eigenstates 
on star graphs, when b is fixed, that are strongly scarred 
by certain short periodic orbits. (Note that on such 
graphs all orbits are unstable.) Our construction exploits 
the properties of the spectral determinant . The spec- 
tral determinant has poles at the points 



b b f 

i=l i=l ^ 



7r/2 -|- im 



Since the derivative of Z(k) is everywhere greater than 
zero, there is exactly one root of Z{k) — between every 
two consecutive poles. 

Given a small e > 0, which will control the quality 
of the scarred eigenstate, wc can find a pole pi in the 
set Pi satisfying the following properties: (a) there is 
a pole p2 from P2 within a distance e of pi and (b) pi 



is approximately equidistant from the two nearest poles 
from Pi, for each i > 2. Due to the ergodic properties of 
the sequence P (assuming that the bond lengths Li are 
incommensurate), the above situation occurs with non- 
zero frequency along the fc-axis. 

Denote the root squeezed between pi and p2 by k' . 
Then cos k' Li is of the order of e when i = 1,2 and is of 
order 1 otherwise. Going back to the eigenstate formula 
©J we see that 



Al,2(fcO 

Mk') 



0{e- 



for i > 2, 



that is, the amplitude of the fc'-eigenstate on the bonds 
1 and 2 is e^^ times stronger than on any other bond. 
By selecting suitably small e one can find eigenstates lo- 
calized on any two given bonds to any precision. Under- 
standably, higher precision leads to a smaller frequency 
of the scarred eigenstates. In fact, the frequency is pro- 
portional to e. 

Since Z{k') = it follows that Ai{k') w A2{k') which 
provides an explanation for the visible singularity at i? = 
1/2 in the difference between P{R) for finite b and its 
limiting form (see figure EJ- This singularity corresponds 
to the eigenstates localized on bonds e and e' such that 
e is picked out by the observable B and e' is not. 

The above construction can be generalized to produce 
eigenstates localised on any number j > 2 of bonds. How- 
ever, once j > 2, the amplitudes on the j bonds are gen- 
erally not equal, which explains the lack of singularities 
at rational fractions other than 1/2. Finally, the singu- 
larities at i? = and 1 correspond to the cases when the 
eigenstates are localized fully outside {R = 0) or inside 
(i? = 1) the V bonds picked out by B. 

The preceding calculations can be made rigorous. We 
defer the details to |2^ . 

In it was suggested that the squares of the co- 
efficients, Ci, of the eigenfunctions of Seba billiards ex- 
pressed in the basis of states of the unperturbed billiard 



(0)\ 



(9) 



are distributed in the same way as the square of the max- 
imum norm on a single bond of a star graph in the limit 
as w ^ 00. This conjecture was supported by numer- 
ical evidence. We extend this analogy to interpret the 
above results in terms of the Seba billiard. Since the 
quantity in Q is similar to a sum of norms of eigenfunc- 
tions on a fraction of bonds, we conjecture that the sum 
of the squares of a fraction a^^ of the coefficients has 
probability distribution P{R)- To elucidate this idea, 
consider preparing a Seba-type system in a randomly- 
chosen eigenstate. The perturbation is then removed in- 
stantaneously, and a measurement of the energy is made. 
What is the distribution (with respect to the choice of 



initial state) of the probability that the measured en- 
ergy is one of a given fraction of the energy levels 
of the unperturbed system? The answer is the distribu- 
tion function in ©. If the eigenfunctions of the billiard 
were asymptotically equidistributed then this probability 
distribution would be a unit step function at R= 1/a. 

Energy levels of a Seba billiard interlace with energy 
levels of the original unperturbed system in much the 
same way that momenta of star graphs interlace with 
poles of the function Z{k). We consider a Neumann bil- 
liard in a rectangle with aspect ratio 7^^^, perturbed by a 
point singularity at the origin. Eigenstates of this system 
can be expanded as 



En 



(10) 



where An is a normalisation constant, the energy levels 
of the Neumann billiard are E^^^ — 47r^7~-^/^(j^ -t- jj"^), 

and IV'i j"*) are the corresponding eigenfunctions. It is 
well known that these unperturbed eigenfunctions are 
localised in momentum space. We therefore expect to 
find states of the Seba billiard that exhibit structures 
analogous to scars in momentum space when their en- 
ergy is between two closely spaced levels of the unper- 
turbed billiard. In fact such states will scar in two di- 
rections in momentum space, corresponding to the two 
unperturbed eigenstates closest in energy to the state in 
question. These scars are supported by families of orbits 
corresponding to tori in the unperturbed system. Note 
however that torus quantisation itself does not apply. It 
is in this sense that the structures are analogous to scars. 

Figure shows the 55"^ state of the Seba billiard 
described above, with the scatterer placed at the ori- 
gin. Although there is no clear localisation evident in 
position representation, the momentum representation 
clearly shows localisation in two directions. 
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